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Abstract. - We consider a mesoscopic normal metal, where the spin degeneracy is lifted by 
a ferromagnetic exchange field or Zeeman splitting, coupled to two superconducting reservoirs. 
As a function of the exchange field or the distance between the reservoirs, the supercurrent 
through this device oscillates with an exponentially decreasing envelope. This phenomenon is 
similar to the tuning of a supercurrent by a non-equilibrium quasiparticle distribution between 
two voltage-biased reservoirs. We propose a device combining the exchange field and non- 
equilibrium effects, which allows us to observe a range of novel phenomena. For instance, 
part of the field-suppressed supercurrent can be recovered by a voltage between the additional 
probes. 



Externally controlled weak links in mesoscopic superconducting circuits have been at the 
focus of interest in recent years Q. The possibility to control the quasiparticle distribution 
by external voltage probes allows tuning the supercurrent through the device (mesoscopic 
SNS transistors). It has been predicted that devices with tunnel junctions and systems 
with good metallic contacts || can enter a peculiar mesoscopic non-equilibrium state at low 
temperatures, which even allows reversing the supercurrent, turning the system into a n- 
junction. This phenomenon has been verified experimentally Jl]]. 

Another phenomenon of high interest in superconducting mesoscopics is the combination of 
ferromagnetic (F) elements with superconductors (S) A strong exchange interaction h in 

the fcrromagnet is expected to suppress the superconducting proximity effect, and hence also 
the supercurrent. (Several recent experiments |p|-pT[ do not confirm this expectation, a fact 
which, at this stage, is not understood.) For weak fields, the supercurrent through a SFS weak 
link and the transition temperature of a SF multi-layer are predicted to oscillate [|l2|-|l4|] as a 
function of the field, or of the width d of the ferromagnet. The latter defines a characteristic 
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energy scale, the Thouless energy, which in the diffusive limit is E^ = hD/d 2 , proportional 
to the diffusion constant D. 

In this article, we show that the non-equilibrium-controlled supercurrent in mesoscopic SNS 
transistors || and the supercurrent in SFS weak links are formally equivalent, although one 
is tuned by varying the distribution function, while the other is controlled by modifications of 
equilibrium spectral functions. Combining the two phenomena, we can recover by an applied 
voltage part of the supercurrent which is suppressed by the exchange field. Thereby, one can 
measure the exchange field in the weak link. 

For definiteness we consider a quasi-one-dimensional system depicted schematically in fig. ^ 
assuming a three-dimensional system with structural changes only in one direction. The 
magnetism in the weak link, or the Zeeman splitting, is accounted for by the energy ah 
of a spatially homogeneous exchange field coupling to the electron spin a = ±1. In the 
diffusive limit, the system can be described by the Usadel equation for quasiclassical Green's 
functions f||7], E^. While the equilibrium results of the present work can also be obtained 
in the imaginary-time Matsubara formalism, we have chosen to use the real-time Keldysh 
technique in order to include also non-equilibrium processes. Then we have 

Ddl6 = - 2i(E - ah) sinh 6 + 2A cosh 6 + —{d x xf sinh 26 (1) 

JE (E, h) = sinh 2 9d xX , d x]E = , (2) 

where 9{E, h, x) and x(-E, h, x) are complex variables parametrising the quasiclassical diagonal 
and off-diagonal Green's function G(E, h, x) = cosh(0) and F(E, h, x) = sinh(0) exp(ix). For 
a system of length d, eq. ([!]) introduces a natural energy scale Et = D/d 2 . Hence, one way to 
tune the relevant energies is by varying the length d. Deep in the superconducting electrodes 
the exchange field or Zeeman splitting vanishes. For simplicity, we assume a bulk BCS solution 
up to the interfaces, 9s = arctan(A/J5), xs — in the superconducting electrodes with 

amplitudes A and phase difference <f> of the order parameters of the two superconductors. 
Furthermore, we assume clean interfaces, and neglect the reduction of Andreev reflection 
expected in spin-polarised systems p^-[l^| . We expect the error due to these approximations 
to be only quantitative (for the latter point, see the discussions below). 

The imaginary part of the conserved 'spectral supercurrent', je, in eq. (||) enters into the 
observable supercurrent as 

jsW = tE»nJ dE(l - 2f(E))Jm{j E (E, ah)}. (3) 

<T=±1 J -°° 

Here f(E) is the distribution function of quasiparticles in the weak link, which in the ab- 
sence of applied voltages reduces to the equilibrium Fermi distribution / cq . Furthermore, 
<7n<t = 2e NaaDcr is the normal-state conductivity for spin a, and Nq^ is the corresponding 
normal-state density of states at the Fermi level. Our approach (eqs. (Q), (||)) assumes spin- 
independent densities of states and diffusion constants. It is valid at low fields h, when the 
variation in the densities of states is small, iVof — Nqi <C (iVo-f + iVoj.)/2. In this case we 
may put gwr= gNi = 9n- The distribution function / in general is obtained from kinetic 
equations jl5| , but for the moment, we assume thermal equilibrium. 

It is instructive to see how the spectral supercurrent Im{j£} depends on energy E and 
exchange fields h. It is plotted in fig. || for a phase difference (f> — it/2 between the supercon- 
ducting electrodes. For h = 0, the function Im{j£} is antisymmetric around the Fermi surface. 
At low energies E < Et, it vanishes until some phase-dependent E c {<f). At larger energies it 
increases sharply, and then decreases exponentially, oscillating between positive and negative 
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Fig. 1 - Schematic picture of the studied S-F-S structure. 

Fig. 2 - Spectral supercurrent for different exchange fields h at <f> — n/2 as a function of energy. 
The exchange fields are expressed in the units of the Thouless energy Et- The variation in the peak 
heights is due to a finite magnitude of the order parameter A — 50Et- 



values. The exchange field shifts the position of the symmetry point from E — to E = ah 
and for a superconducting gap A of the order of h, distorts the symmetry. Since A serves 
as an upper cutoff, which is not shifted, the overall magnitude of the spectral supercurrent 
decreases when h becomes comparable to A. 

In equilibrium we have 1 — 2f(E) = tanh(E/2T). This term and the sum of the spectral 
supercurrent J2 a jE(E, ah) are antisymmetric around E = 0. Hence, for the discussion of the 
total supercurrent jg we can concentrate on the part E > 0. At low T < Et the supercurrent 
js is given by an alternating sum over the decreasing areas under the oscillating function 
Imjjs} measured from the i?-axis (see fig. ||). At h = 0, the positive first term dominates the 
sum and yields a large supercurrent js. Increasing h shifts the negative peak from E < to 
E > 0, hence decreasing js, and even reversing its sign. At finite temperature, the low-energy 
part, E < T, is effectively cut off, hence js decreases in amplitude. This result is illustrated in 
Figs. H and ||, where js {(j) = 7t/2) is plotted as a function of different exchange fields at different 
temperatures and for different bulk order parameters A. Analogous results can be obtained 
for a constant exchange field by varying the distance d of the superconducting reservoirs and 
through it the Thouless energy Et- 

In the regime where js (</> = 7r/2) is negative, the junction forms a so-called l 7r-junction' fl~2|| , 
since the ground state of the system, with no supercurrent flowing between the two super- 
conductors, is reached for a phase difference equal to n. The supercurrent-phase relation for 
different exchange fields is plotted in fig. [5], showing the crossover from an ordinary behaviour 
to a 7r-state. A closer analysis of fig. || shows that the precise value oih/ET where the crossover 
occurs depends weakly on temperature, since higher values of T smoothen the oscillations be- 
tween positive and negative contributions to js. At h = 0, js(</> = t/2) is positive for any T. 
It remains positive as long as the thermal energy dominates over the exchange, T S> h. With 
increasing T the cross-over to a 7r— junction is shifted towards higher fields. This dependence 
was probably observed in ref. jOI. It is an alternative way to verify the current reversal to 



4 



EUROPHYSICS LETTERS 




10 



20 
h/E n 



30 



40 



10 



20 
h/E n 



30 



40 



Fig. 3 



Fig. 4 



Fig. 3 - Supercurrent js(4> = 7r/2) as a function of exchange field h/Ei through the structure depicted 
in fig. [j] for different temperatures T/Et- The superconducting order parameter A = 1000-Et- 

Fig. 4 - SFS supercurrent js{<t> = tt /2) as a function of exchange field }i/Et for different values A of 
the superconducting order parameter at T — 0. 



what has been discussed in previous proposals, where typically one requires many different 
samples with varying widths Q but otherwise equal parameters. The crossover is illustrated 
in fig. H for a few values of h/Er- 

By shifting the variable of integration E in eq. (J^j by ah one finds 

k{4>) = E r dE ^ - + "h))Tm{j E (E)}. (4) 

The shifted distribution function / = 1/2 ^ CT f cq (E + ah) has the same form as the two-step 
distribution function measured in ref. pO| j. There it appeared as the solution of a kinetic 
equation in the centre of a diffusive metal between two normal probes with voltage eV = ±2h 
in the limit where the inelastic scattering length is longer than the distance between the two 
normal reservoirs. The spectral supercurrent in general still depends on the exchange field 
via the boundary conditions. However, if the superconducting gap A is much larger than the 
exchange field, A 3> h, this dependence can be neglected. In this limit, the supercurrent in 
the presence of an exchange field is the same as for a non-equilibrium distribution four-probe 
structure described in ref. || (see fig. [7] for a schematic picture) . 

It is interesting to note that this behaviour of the diffusive-limit supercurrent as functions 
of the exchange field and the external potential is very similar to the supercurrent through 
a multi-probe structure in the clean limit. This limit has been described by Dobrosavljevic- 
Grujic et at. 0| for the ferromagnetic two-probe case and by van Wees et at. p2| including a 
voltage in a non-magnetic three-probe setup. In this case, the supercurrent is carried by the 
Andreev levels, whose energies are controlled by the exchange field Et] , and whose occupation 
can be tuned by the voltage [^] . With both parameters, for example, the system can be driven 
into a 7r-state. 

We can combine the effects of exchange field and non-equilibrium distribution ||, 23 1 by 
considering the structure in fig. fj], where the magnetic material is placed between supercon- 
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Fig. 5 - Supercurrent js(0) as a function of phase for different exchange fields 1%/Et in the regime 
where the crossover from the ordinary 0-state to the 7r-state occurs for the first time. Here, T = 
and A = IOOCLEt. 



Fig. 6 - Crossover from the 7r-state (js(4> = tt/2) < 0) to the 0-state (js( 
of temperature T for a few values of the exchange field }i/Et- 



7r/2) > 0) as a function 



ductors and normal voltage leads. Here, the distribution function is 

f(E, y) = {\- |)r*(£ + eV/2) + (i + j-)f eq (E eV/2), (5) 

exhibiting the two-step form observed by Pothier et al. p0[ . In this case, if the superconducting 
reservoirs are located around y — and provided A>/i, eV, the observable supercurrent can 
be written as a sum of four terms, 

JsW = !°° d£(l - 2n(E))lm{j E (E -h- eV/2) + j E (E -h + eV/2) + 

8 J-oo (6) 

je(E + h- eV/2) + JE (E + h + eV/2)}. 

For example, if the potential is exactly twice the exchange field, eV = 2h, due to the anti- 
symmetry of Imjj^;}, we have 

M<t>) = \ (j! FS (<f>,o) + ji FS (^h)) « i Js SFS (0,o). (7) 

The latter approximate equality holds if h ^> Et. Here, j| F (</>, /i) is the supercurrent through 
the SFS structure with the exchange field h in the weak link. Hence, one can use the external 
potential to recover half of the zero-field supercurrent. This is illustrated in fig. |[ where 
the supercurrent of the four-probe structure is plotted as a function of voltage with a few 
magnitudes of fields. 

The results summarised by eq. (^) provide a way to measure the exchange field and at the 
same time to explore the applicability of the simplified model for the ferromagnet used here 
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Fig. 7 Fig. 8 



Fig. 7 - Four-probe setup for studying the non-equilibrium effects on the supercurrent. It is assumed 
that L 3> ds and that the superconductors lie in the middle of the normal wire (y — 0) so that the 
distribution function has the two-step form between the superconductors. Furthermore, we expect 
that the four-probe setup does not notably alter the spectral supercurrent obtained from a quasi-one- 
dimensional calculation (for a more detailed discussion, see refs. 

Fig. 8 - Supercurrent js(4> = tt/2) of the four-probe structure at different fields as a function of the 
voltage V between the normal probes. In the calculations for the main picture, the magnitude of the 
order parameter was set to A = 100i?T, and at the inset, A = 10-Et, thereby showing that even when 
A is of the order of h and eV, a local maximum is obtained at eV = 2h. 



and previously When the voltage-dependent supercurrent js(V) reaches a maximum, 

eV should equal 2h. Deviations could occur as, for instance, this model neglects the band- 
structure effects [pij important in the ferromagnets. Also, to be able to measure the actual 
supercurrent through a typical ferromagnet with Curie temperature Tq u ^> A, the ratio h/Ex 
has to be made small by fabricating very thin weak links. Moreover, our assumption of the 
diffusive regime requires d 3> l c \, and a quantitative agreement cannot be expected for thin 
structures. Finally, due to the strong electron-electron interactions in ferromagnets, producing 
a short inelastic relaxation length, the normal probes should be fabricated rather close to each 
other to obtain the two-step form for the distribution function. For conventional ferromagnets 
the exchange field is large, which would correspond to enormous voltages. However, we expect 
that our model is approximately valid for setups constructed from ferromagnetic alloys with 
Tq u of the order of the superconducting critical temperature fl9[| , or in situations where h can 
be related to the Zeeman splitting in magnetic fields much weaker than the superconducting 
critical field. 

In summary, we have calculated the supercurrent through a ferromagnetic weak link as a 
function of the exchange field in the ferromagnet. In the calculations, the Keldysh technique 
was used to provide a description of non-equilibrium effects. We found that when A 3> h, the 
problem is formally equivalent to the four-probe measurement of the supercurrent through a 
normal-metal weak link. Furthermore, we showed that applying a non-equilibrium potential 
in the transverse direction, one can recover half of the supercurrent of a ferromagnet with an 
exchange field h ^> _Ex, as compared to the supercurrent in the absence of h. 
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